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One of the most fundamental problems in quantum many-body systems is the identification of 
a mean field in spontaneous symmetry breaking which is usually made in a heuristic manner. We 
propose a systematic method of finding a mean field based on the Lie algebra and the dynamical 
symmetry by introducing a class of symmetry broken phases which we call /r-symmetry breaking. 
We show that for /r-symmetry breaking the quadratic part of an effective Lagrangian of Nambu- 
Goldstone modes can be block-diagonalized and that homotopy groups of topological excitations 
can be calculated systematically. 


I. INTRODUCTION 


Spontaneous symmetry breaking (SSB) has long 
played a pivotal role in our understanding of Nature [I|. 
Examples include ferromagnetism superconductivity 
Q, Bose-Einstein condensation chiral symmetry 

breaking @,0 , and unification of the fundamental forces 
Q. Both static and dynamic properties of a symme¬ 
try broken phase can be described by the corresponding 
mean field, which is usually found in a heuristic manner. 
The identification of the mean field amounts to that of 
an order parameter or that of an operator that supports 
a long-range order (LRO) in quantum field theory [9|. 

In this paper, we propose a systematic method of find¬ 
ing mean fields of quantum many-body systems based on 
the Lie algebra and the dynamical symmetry. The dy¬ 
namical symmetry has achieved a remarkable success in 
few-body systems for finding the atomic spectrum of hy¬ 
drogen [l^ [ll| and collective excitation spectra of nuclei 
[l^ . Here we apply the dynamical symmetry to a partic¬ 
ular class of broken symmetry systems in which the mean 
fields are described in terms of the weight vector in the 
representation of the Lie algebra. Since the weight of the 
Lie algebra is often labeled by the Greek letter fi, we refer 
to such symmetry breaking as /i-symmetry breaking. We 
show that for /r-symmetry breaking the quadratic part 
of an effective Lagrangian of Nambu-Goldstone modes 
(NG modes) can be block-diagonalized and that homo¬ 
topy groups of topological excitations can be calculated 
systematically. By applying this method to a U{N)- 
symmetric system which has recently been realized in an 
ultracold atomic gas EH, we show that a large class 
of symmetry broken phases can be described in terms of 
/i-symmetry breaking. 

This paper is organized as follows. In Sec. m we in¬ 
troduce the concept of /i-symmetry breaking and derive 
mean fields by combining it with the dynamical sym¬ 
metry. In Sec. uni mean fields of /i-SB are derived 
through minimization of energy functionals constructed 
from the underlying Lie algebra. In Sec. IIVI we show 
that the quadratic part of an effective Lagrangian of 
Nambu-Goldstone modes can be block-diagonalized for 


/i-symmetry breaking. In Sec. El we show how to system¬ 
atically calculate homotopy groups of topological excita¬ 
tions for /i-symmetry breaking. In Sec. ED we apply our 
method to a I/(iV)-symmetric system. In Sec. IVIIl the 
cases of higher-dimensional representations are discussed 
from the standpoint of /i-symmetry breaking by using 
examples of spin-2 Bose-Einstein condensates (BECs) 
[lB-Ii3 s-iid spin-1 color superconductors [3, m. In 
Sec. IVIIIi we conclude this paper. In Appendix „ we 
prove some formulas on homotopy groups used in Sec. El 


II. /i-SYMMETRY BREAKING 

We consider a quantum field theory whose symmetry 
group G is described by a finite-dimensional unitary rep¬ 
resentation R: 

( 1 ) 

ij 

where {4>i}i is a set of fields of particles, is an ele¬ 
ment of the Lie algebra g = {Ta}a=i constituted from 
finite-dimensional Hermitian matrices of the representa¬ 
tion R, respectively. Here, d is the dimension of G and 
taS are real parameters. We denote the Noether charge 
associated with the generator Ta by Qt„ ■ For the present 
discussion, we do not need to specify quantum statistics 
of particles and the system can be defined either on a 
lattice or in continuous space. We assume that the Lie 
algebra g of the symmetry group is the direct product of 
a simple compact Lie algebra g and u(l) = {xl\x G R}: 

g = g©u(l), (2) 

where / is the identity matrix. The particle-number oper¬ 
ator is the Noether charge associated with the generator 

I. 

The key ingredient in the following analysis is the 
quadratic Gasimir invariant defined by 

Cl{v)-.= J2{v\Ta\v)\ (3) 

a=l 
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where |v) and d = d—1 are a vector in the representation 
R and the dimension of 0 , respectively. Let be 

the Cartan subalgebra of g, i.e. the maximal commuta¬ 
tive subalgebra of g, where r = rank g is the rank of the 
Cartan subalgebra. Let |/x) be a weight vector which is 
a simultaneous eigenstate of 

iLb|/x) =^b|/x) (6 = ,f), (4) 

= ,Hf), (5) 

where t denotes the transpose. The highest weight Ifin) 
is the weight vector that maximizes the expectation value 
of the quadratic Casimir invariant Cl: 


As is shown in Eq. (l32l) . Eq. © implies that the sum of 
the weight vectors within each unit cell vanishes. 

The derivations of Eqs. ([5]) - (ITT]) will be shown in the 
following section. We call these four types of symmetry 
breaking /i-symmetry breaking (/r-SB), which is charac¬ 
terized by the combination of the highest or zero weight 
and ODLRO or DLRO. Prototypical examples of /x-SB 
include the ferromagnetic phase and the polar (antiferro¬ 
magnetic) phase of a spin-1 BEC [2l|,[22, classical ferro- 
magnets (FMs), and classical antiferromagnets (AEMs) 
(see Table m. 

The order parameter of a spin-1 BEC is described by 
a three-dimensional complex vector 


d 

a—1 

d 


( (</'!) \ 

{4>) = (^o) , (13) 

V (</>-!) ) 


= max 




( 6 ) and the symmetry group G is C(1) x S'0(3). The Cartan 
generator of SO{3) is the S'z-operator defined by 


The zero weight |/xo) is the weight vector that has the 
zero eigenvalue of Eq. o and therefore minimizes the 
expectation value of C® to 0 : 


/I 0 0 \ 

= p 0 0 . (14) 

VO 0 -ij 


C|(/xo) = 0. 


(7) The quadratic Casimir invariant of this system is 


An irreducible representation of g is uniquely determined 
by its highest weight fin [ 1 ^, and we denote a complete 
set of weight vectors belonging to fin as Wlfin]- 

We now introduce the concept of /i-symmetry break¬ 
ing. Let {4>) be an order parameter. If (0) trans¬ 
forms in a low-dimensional representation of the sym¬ 
metry group G in the presence of off-diagonal long-range 
order (ODLRO), {4>) will be shown to take either of the 
following forms: 

(</>) = Ifin), ( 8 ) 

{cj}) = l/io). (9) 

If the lattice of space, which we denote by L, is free from 
frustration in the presence of diagonal long-range order 
(DLRO), the mean-field ground state IG^) will be shown 
to take either of the following forms: 

|G5)=(g)|/iP, (10) 

i£L 

|G5) = (g) (g) l/i^ , (11) 

u^lA iGu 


where |/ii)j, u, and U denote a simultaneous eigenstate 
of at lattice site i, a unit cell of an ordered state, 

and the lattice constituted from the entire set of the unit 
cells, respectively. In Eq. m, the set {Hi}i is chosen 
so that the expectation value of G| within each unit cell 
vanishes (see Eq. (l32l) and the following explanation for 
detail): 


(C!) = E 




.iGu 


G|‘’^^^ = 5-5, (15) 

where S = {Sx, Sy, Sz) is the vector of spin operators in 
the Cartesian representation. In the ferromagnetic phase, 
the order parameter has the form 


{4>)fm = 



(16) 


This is the eigenstate of Sz with eigenvalue 1 and maxi¬ 
mizes the expectation value of = S ■ S: 

(G|“('V = {S) ■ {S) = 1. (17) 


Thus, the ferromagnetic phase of a spin-1 BEC is char¬ 
acterized by /i-SB with ODLRO and the highest weight 
In the polar phase, the order parameter has the 

form 


(^)polar 



(18) 


This is the eigenstate of Sz with eigenvalue 0 and mini¬ 
mizes the expectation value of G|°^^V 

(G|“^^V = 0. (19) 

Thus, the polar phase of a spin-1 BEC is characterized 
by /x-SB with ODLRO and the zero weight /xq. 

For both the classical EM and the classical AFM, the 
quadratic Casimir invariant is again the square of a spin 


= 0. (12) 


( 20 ) 
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TABLE I. Examples of four types of /r-symmetry broken 
phases. Here, ij,h and fig represent the highest-weight and 
zero-weight vectors, respectively. The order parameters for 
ODLRO are given by Eqs. (O and m and the ground states 
for DLRO are given by Eqs. ra and dj. 



ODLRO 

DLRO 

/Xif 

ferromagnetic spin-1 BEC 

ferromagnet 

/Xo 

polar spin-1 BEC 

antiferromagnet 


Let S be the spin quantum number. The highest- 
weight state IfJ-n) corresponds to the eigenstate of 
with the highest magnetic quantnm number ruz, i.e., 
iH/f) = {mz = S). The mean-field ground state |G5 ')fm 
of the classical FM represents a uniform alignment of the 
highest magnetic quantum-number state \mz = S) on 
every site of the lattice L: 

|G5)fm = (g) |m, = 5),. (21) 

i^L 

Therefore, the classical FM is characterized by /r-SB with 
DLRO and fin- On the other hand, the mean-field 
ground state |GS')afm of the classical AFM is described 
by a uniform alignment of the highest magnetic quantum- 
number state \mz = S) on the sites of one sublattice La 
and that of the lowest magnetic quantum-number state 
\mz = — S') on the other sublattice Lb- 

|GS)afm = (g) {\mz = S)uA O \mz = -S)ub) , (22) 

where ua and ub indicate the sites in each unit cell u 
belonging to La and L^, respectively. Within each unit 
cell, the total magnetization vanishes. Therefore, this 
phase is characterized by /x-SB with DLRO and /xq- 


III. DERIVATIONS OF MEAN FIELDS OF fi-SB 

In this section, we show that mean fields described by 
the highest or zero-weight states are obtained through the 
minimization of an energy functional constructed from 
Casimir invariants of the underlying Lie algebra g. 

We first discuss the case of ODLRO. In this case, the 
energy functional can be obtained in a manner similar to 
the case of the dynamical symmetry [ni[il,|23. How¬ 
ever, as shown later, for ^-SB we can block-diagonalize 
the quadratic part of an effective Lagrangian of NO 
modes and systematically calculate homotopy groups of 
topological excitations, neither of which can be done 
from the dynamical symmetry alone. For the case of 
the lowest-dimensional representation such as a BEG in 
degenerate TV-component bosons the mean-field energy 
functional V{{4>)) can be expressed up to the fourth or¬ 
der in {cj)) as 


where c and cq are real constants. Since any minimizer 
(</>) of V{{cf>)) in the lowest-dimensional representation 
can be transformed into \hh) by an appropriate element 
of G, the energy functional is minimized for 

(</») = IfJ-H) ■ (24) 

For the case of the next lowest-dimensional represen¬ 
tation such as a spin-1 BEG [2l|, and an s-wave 
superfluid in degenerate V-component fermions [23 - [^ 
V{{4>)) can be expressed in terms of |(</>)P and the 
quadratic Gasimir invariant as 

E(( 0 )) = -cm\^+coim^+c^cim), (25) 

where c, cq, and ci are real constants. For the case of 
a ferromagnetic interaction between condensate particles 
with Cl < 0 , the energy functional is minimized for 

( 0 ) = M ■ (26) 

For the case of an antiferromagnetic interaction with Ci > 
0 , the energy functional is minimized for 

(</>) = IHo) , (27) 

if the representation includes the zero-weight state |/xo). 
This condition is satisfied for low-dimensional Lie al¬ 
gebras such as su(2) and so(3). The general case of 
0 = su(TV) will be discussed in Sec. |VT1 Since the en¬ 
ergy functional can be written in neither form of Eq. (EH) 
nor Eq. ()25|) in higher-dimensional representations, the 
ground states are no longer described by ^-SB. Such a 
case will be discussed in Sec. IVIII 

We next discuss the case of DLRO. We assume that 
particles are placed on a lattice L whose geometry is free 
from frustration. A prototypical example of DLRO is 
described by the Heisenberg-type Hamiltonian 

d 

H=-jJ2Y.^aATa,j, (28) 

{ij) “=1 

where {i,j) represents a pair of nearest-neighbor sites i 
and j, and {TaA}i=i is a set of generators of a simple 
compact Lie algebra g on site i. For the case of a ferro¬ 
magnetic interaction with J > 0 , the ground state [GS") 
is written as 

|G5) = (g) l/xxx),, (29) 

iei 

which coincides with Eq. (nni). For the case of an an¬ 
tiferromagnetic interaction with J < 0 , it follows from 
the frustration-free assumption that the mean-field clas¬ 
sical ground state \GS) is obtained by a tensor product of 
the state on site i that minimizes the interaction energy 
with its neighboring sites [2^. Therefore, the mean-field 
ground state can be written in terms of a site-factorized 
wave function as 

|G^) = (g)(g)lH4, 

u£U i^u 


vm) = -om\^+c,m\\ 


(23) 


(30) 
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where u, and lA denote a simultaneous eigenstate 

of {Hi,Yb=i lattice site *, a unit cell of an ordered state, 
and the lattice constituted from the entire set of the unit 
cells, respectively. In Eq. (1501) . two weight vectors and 
^Jbj at nearest-neighbor sites i and j should satisfy 

(/Xi|/Xj) = min{(/Xfc|/Xi) \^lk, £ WIiih]} (31) 


to minimize the nearest-neighbor interaction energy 
Sa=i Ta,iTaj. We note that we do not impose Eq. (fT^ 
on Eq. (pl) so far. Since the expectation value 
(fii \Ta\Hi) vanishes for any off-diagonal matrix Ta, 
which is nothing but the raising or lowering operator of 
the Cartan canonical form (see Sec. IIV Al) . the expecta¬ 
tion value of C® within each unit cell u can be calculated 
from Eqs. @ and (1301) by considering the contribution 
from the Cartan generators {Hb}Yi alone. We thus ob¬ 
tain 


b—1 Li£u 


(Cl) = E 

f 

= E 


^ ^ 

i^u 


i^u 


(32) 


where ||a;|| denotes the magnitude of the vector x. If 
there exist a pair of weight vectors with opposite direc¬ 
tions such as the case of SU{2) (see Fig. [T](a)), Eq. (1^ 
is satisfied and the sum of the weight vector in Eq. (15^ 
vanishes within each unit cell. Thus, the ground-state 
wave function satisfies Eq. m- However, in a larger 
group such as 51/(3), it is known that there do not ex¬ 
ist two weight vectors with opposite directions and that 
there is more than one pair of weight vectors that satisfy 
Eq. (15X1) (see Fig. [I](b)) [2^. For such cases, the mean- 
field ground states represented by Eq. dsni) are degenerate 
[^ . Therefore, we have to consider higher-order contri¬ 
butions arising from quantum fluctuations to determine 
the ground state. This can be done by using the flavor- 
wave theory [29l - l3^ . The Hamiltonian iXfw of quantum 
fluctuations is given as follows [3l|, : 

Hi-w = E! (^^) 

(iJ) 


where 



FIG. 1. (a) Weight vectors of the two-dimensional represen¬ 
tation of the su(2)-Lie algebra. The weight vectors are two 
one-dimensional vectors izi = l /%/2 and t '2 = — l/v^, which 
have opposite directions with the same magnitude l/%/2. In 
this representation, the pair of the weight vectors that satisfy 
Eq. (I31II is uniquely determined to be (i/i,r' 2 ). (b) Weight 
vectors of the three-dimensional representation of the su(3)- 
Lie algebra. The weight vectors are three two-dimensional 
vectors ui = (l/%/ 2 , l/p), 1/2 = (—l/v^, l/v/6), and 1/3 = 
(0, —■\/2/3) from the origin to each of the three apexes of the 
equilateral triangle. The sum of these three weight vectors 
vanishes. Three pairs of the weight vectors, (jzi, 1 / 2 ), {i' 2 , 1 ^ 3 ), 
and ( 1 Z 3 , u\), satisfy Eq. (13111 . In both (a) and (b), the weight 
vectors are normalized so as to satisfy Eq. (IIO 8 II . 


follows from Eq. (1551) that 

= (36) 

and hence we have 

+ Mfc (37) 

for any pair of sites i and k that have a common nearest- 
neighbor site j. Thus, the ground state must satisfy both 
Eqs. (IHTl) and (1571) . Equation (1571) implies that two weight 
vectors, fn /Xj, must give the minimum value of the 
inner product (/Xi|/Xj) for any pair of nearest-neighbor 
sites i and j. Equation (EZl) implies that two weight vec¬ 
tors, Hi s-nd Hk, must be different for any pair of sites i 
and k that share a common nearest-neighbor site. As a 
consequence, the sum of the weight vectors within each 
unit cell in Eq. (I55|) tend to vanish and an ordered state 
with the vanishing expectation value of G| within each 
unit cell is favored. Thus, the ordered state is character¬ 
ized by ^-SB with ho- 


IV. EFFECTIVE LAGRANGIAN OF NG 
MODES FOR /x-SB 


Aij — + b^nji, (34) 

and is the annihilation operator of a boson with 

with weight vector Hi at site j. The expectation value of 
i7fw is minimized when 

(G5|A^Ay |G5) = 0 Ay |G5) = 0 for Vx, j. (35) 

Let i,j and j, k be pairs of nearest-neighbor sites. Then, 
i and k have the common nearest-neighbor site j. It 


A. Cartan canonical form and the generalized 
magnetization 

The analyses of NG modes and topological excitations 
in ^-SB can be done conveniently in terms of a special 
basis of the Lie algebra g known as the Cartan canonical 
form [ 23 : 




( 38 ) 





















where R and / indicate the real and imaginary parts of 
the raising operators E±c<l of the Cartan canonical form: 
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:= ^ 


E_ 




EL := 


Err — E—r 


V2i 


(39) 


The Cartan canonical form (l38l) is a generalization of 
the basis of the su(2) Lie algebra 5”^}}, which 

decomposes the generators into the diagonal matrices 
{-^b}&=i fhe off-diagonal ones {i?^, i?^}aG-R+i where 
a. is an f-dimensional real vector known as a positive root 
and i?+ denotes the set of all the positive roots. Here the 
positive roots distinguish different su(2)-subalgebras of g. 
Defining Ha for a S R’’ by 

f 

Ha := Yr ^bHb, (40) 

b=l 


where 


a — *(«!, a2, • • ■ , cif), 

we obtain a triad Sa for a positive root a as 

Sa^{E!^,Ei,Ha). 


(41) 


(42) 


The triad Sa forms an su(2) subalgebra of q analogous 
to the su(2) spin algebra. In fact, it satisfies the following 
commutation relations: 


[E^,Ei] 

= i{a, a)Ha for Va G i?+, 

(43) 

[ELHa] 

= i{a, a)Ea for Va G i?+, 

(44) 

[Ho.,E^] 

= i{a,a)Ea for Va G R+, 

(45) 

[Ha, Ha' 

] = 0 for ya,ya' G R+- 

(46) 


We refer to Sa as a generalized magnetization in analogy 
with the su(2) spin S. As shown later, the textures of 
NG modes and topological excitations are described in 
terms of the generalized magnetization Sa- 


B. Broken generators for /r-SB 

To determine the quadratic part of an effective La- 
grangian of NG modes, we first prove the following the¬ 
orem on the spaces of broken generators g/f), where (] is 
the Lie algebra of the remaining (unbroken) symmetry of 
the state. 

Theorem 1 Consider /r-SB phases and define the sets of 
positive root vectors Rh, Rq^, and i?Q by 

Rh ■= {a e R+\{a., fin) ^ 0}, (47) 

R^ ■.= \J{aGR+\\fi^-a)GW[fiH]}, (48) 

iGu 

R^^ ■.= {aGR+\\fio-a)&W[fiH]}. (49) 


Then the bases of the space of broken generators g/h are 
given by 


basis (g/()) 
({ES,Ei\ae 
\{ES,EL\ae 
]{ES,Ei\ae 
[{E!^,Ei\ae 


Rh} 

for 

R^} 

for 

Rh] U {/} 

for 

RqOD} U {/} 

for 


DLRO and fin', 
DLRO and /.to; 
ODLRO and fin', 
ODLRO and /xq. 

(50) 


Proof: 

First, we consider the DLRO with /j-h- From the mean- 
field ground state in Eq. m, unitary transformations 
generated by Ea’^ leave the ground state unchanged up 
to a global phase if and only if 

M = 0 

E+a \P-h) = E-a Iftif) = 0. (51) 

Since jj-h is the highest weight, E^a |/^ff) always van¬ 
ishes [^, so that Eq. (ICT) is equivalent to 

E-a Iftff) = 0 \fiH — ol) ^ (52) 

which, in turn, is equivalent to 


a. ^ Rh- (53) 

This equivalence can be shown as follows. Erom Eq. (EH), 
we obtain 


0 = \^JiH) = {(y.a)Ha\^lH) 

= (a,a)(/Xi/,Q:) l/Xi/). (54) 

Since {a, a.) ^ 0 and Ihh) 0, we obtain {fiH,ot) = 0 
and hence Eq. (1551) from the definition of Rh - Conversely, 
we can derive Eq. dMl) by assuming Eq. (l53l) . The Weyl 
reflection of + ot with respect to a is 

, , \ r,iP‘H + OL,a.) 

[HH + a)-2 - - --- ol = hh-ol. (55) 

(a, a) 

Since Ihh) is the highest weight, \fiH + ct) is not a weight 
vector, nor is \fiH — ol)- Thus, Eq. (1551) is obtained. Uni¬ 
tary transformations generated by and / change 

the ground state IGS”) only by a global phase factor. 
Therefore, these generators are not broken ones, which 
completes the proof of the first row of Eq. dsol). 

Second, we consider the DLRO with hq. Erom the 
mean-field ground state in Eq. dm, unitary transforma¬ 
tions generated by Eff'^ leave \GS) invariant if and only 
if 


Ea^ l/Xi) = 0 for Vf S M. (56) 

This condition is equivalent to a ^ From the dis¬ 

cussions similar to the case of DLRO and fin, all of the 
Cartan generators are not broken ones, which completes 
the proof of the second row of Eq. (l50l) . 

Third, we consider the ODLRO with the highest 
weight, i.e. {(p) = {(j-h)- Unitary transformations gener¬ 
ated by EE^ leave \GS) invariant if and only if Eq. (1551) is 
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satisfied. From above discussion, this condition is equiv¬ 
alent to ct ^ Rh- Unitary transformations generated by 
and I change only by a phase factor. This 
phase shift can be eliminated by taking a linear combi¬ 
nation of and I except for the direction of I, 

which completes the proof of the third row of Eq. (|50|. 
The proof of the fourth row can be given similarly by 
replacing fin by /xg and using Eq. ((5^ . Thus, the proof 
of Theorem [T] is completed. 


C. Effective Lagrangian of NG modes 

We now derive the quadratic part of an effective La¬ 
grangian of NG modes. We note that in non-relativistic 
systems the type-2 NG mode with a quadratic dispersion 
is allowed, in contrast to relativistic systems where only 
the type-1 NG mode with a linear dispersion is allowed 
[^ . Let Nb ■= dim{G/H) be the number of the bro¬ 
ken generators. It has been shown that the numbers of 
type-1 and type-2 NG modes, ni and 712 , can be deter¬ 
mined only from the ground state |GS') and the set of the 
Noether charges associated with the broken generators of 
the symmetry group {Ta'}^f,i as follows (the prime in¬ 
dicates that the generators are broken ones) [3, : 

( ni+ 2n2 = dim(G/iL); 

[ 712 = ^rankp, ^ ^ 

p,,b, =-i{GS\[Qt^,,Qt,,]\GS), (58) 

where p = {Pa'h'}^^,^is a Kostant-Kirillov symplectic 
form (K-K form) which is employed in Ref. [s^ 

in the context of NG modes. The basis of the Lie alge¬ 
bra that block-diagonalizes Eq. ([5^ is constituted from 
a set of canonically conjugate pairs and generates type- 
2 NG modes [H, [S^. While the types and the num¬ 
bers of NG modes can be found from Eqs. dSZl) and (I58p . 
the dynamics and the corresponding broken generator of 
the NG mode cannot be determined from them since the 
quadratic part of an effective Lagrangian is not diago¬ 
nalized in Refs, [l^. Here, we show that for p-SB, 
the K-K form and the quadratic part of an effective La¬ 
grangian in the fields of NG modes can be simultaneously 
block-diagonalized in terms of the Cartan canonical form. 
Moreover, NG modes are classified into three categories 
according to their dynamics. In fact, we can prove the 
following theorem. 


Theorem 2 Gonsider a ^-SB in a non-relativistic sys¬ 
tem. Let TT^, TT^ and ttj be the fields of the NG mode 
generated by the broken generators and I and 

define Hq, by Hq, := ir^ + iir^. Then, the quadratic parts 
of effective Lagrangians £eS in Hq, and ttj can be block- 


diagonalized in terms of n„ and tt/ as follows: 

f EaGflH -^pre DLRO and PH\ 

jr ^ I Saefl? for DLRO and /xg; 

I Gpha + -^pre ODLRO and PH-, 

[Upha + ZaGiJon ^osc ^r ODLRO and /xg, 

(59) 

- n^5in„) + 6„|vn„|2, (eo) 

= + (61) 
Gpha = + gi{dtTTiY, (62) 

where Pa,ba,ba, gi, gi are real constants. 

Proof: Gonsider the quadratic part of an effective La¬ 
grangian in the fields of NG modes in a non-relativistic 
system. Then, the most general form of it can be written 
as follows 

Nb . ^ 

Ueff — ^ ^ f Pa'b' '^a' dt^h' ~zga'h' 

a',b' = l ^ 

+ ^5a'b'V7ra' • , (63) 

where p is the K-K form defined in Eq. (l58)) . 
is the set of the fields of NG modes associated with the 
set of the broken generators and fja’b' and ga’v 

are real constants. 

We first block-diagonalize the quadratic part of an ef¬ 
fective Lagrangian. Here we consider the case of ODLRO 
and fiH- The proof for the other cases can be made in 
a similar manner. The quadratic forms constructed from 
7 r/,n^, and Hq are the following six terms: 

(7r7)^,7r7n);,,7r/n„; 

n„n^,n;n*,n„n;3. (64) 

Let be a basis of (f — l)-dimensional subspace 

that is orthogonal to pn- Under the unitary transforma¬ 
tion generated by 

f—1 

+ (65) 

777 is invariant since the generator I commutes with Hs- 
Using the commutation relations of the Gartan canonical 
form 

[E^,H4=z{cx,s)Ei, (66) 

[Ei,E,] =-^(a,s)E^, (67) 

we have 

q-^h,(eR + = e*(^’“)(Uf -h iE^), (68) 

q-^Hs(eR _ = e-*(®’“)(Uf - iE^). (69) 

Therefore, the corresponding fields Hq, and H^ transform 
under the same unitary transformation into Hcte®^®’™^ 
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and respectively. Among the quadratic forms 

in Eq. (IMl) . only and are invariant under the 

transformations generated by Hg for any s G R’’. Thus, 
£eff in Eq. (1551) can be written as 

£eff = PlTTldtTTl + 5/(V7r/)^ + gi{dtTTiY 
ol^Rh 

+6c|vn„p + 6c|atn„|2). (70) 

where pa,p'^,ba, and 5 q, are real constants. The first, 
fourth, and fifth terms correspond to the first term in 
Eq. ((63ll and hence to the K-K form in Eq. ([58|) . Erom 
Eq. (I58ll . we have 

Pi = -i{[Qi,Qi]) = 0, (71) 

Pa = —i{[QE^ + 'I-Qe^^Qe^ ~ 

= -HQhJ, (72) 

P'a = ~ ^Qe^^Qe^ + 

= -Pa- (73) 

Thus, we obtain the block-diagonalized form of the 
quadratic part of an effective Lagrangian in 11^, tt/, : 

£eff = 5/(V7r/)^ -(- giidtTTif 
ol^Rh 

+5„|vn„|2 + 6„|9tn„p], (74) 

= -2i{QHj. (75) 

Here, Rh represents the set of the positive root a vectors 
associated with the broken generators. Eor the ODLRO 
(DLRO) with fiQ, the quadratic part of an effective La¬ 
grangian is obtained by replacing Rh into (Rq)- 
Next, we calculate {QhJ) in Eq. (1751) . For /i-SB with 
MOj {Qho) vanishes for any positive root vector ot be¬ 
cause the expectation value of any generator also van¬ 
ishes. Therefore, the term Pa{^adt^a — n^cltHct) in 
Eq. (l74l) vanishes and hence we obtain the fourth row of 
Eq. (1551) . On the other hand, it follows from the defini¬ 
tion of ^-SB for plh in Eqs. ([8]) and (|9]) that 

{Qho) = Nc{hhtOl) ^ 0 for OL G Rh, (76) 

where Nc is the number of particles that contribute to 
the LRO. The term bcy.\dtTla\^ in Eq. (I74p is much smaller 
than the term PaiJ^adt^^ — n^i9tn„) at low energy be¬ 
cause the former is the second order in dt while the latter 
is the first order in dt- Therefore, the term ba\dtl\-a\^ 
can be ignored and hence we obtain the third row of 
Eq. (l59l) . For the case of DLRO, the terms involving tt/ 
are absent since I is not a broken generator from Theo¬ 
rem [T] Therefore, we obtain the first and second rows of 
Eq. (155)) . which completes the proof of Theorem [5] 


D. Three types of NG modes in /r-SB 

From Theorem[2l we find that three types of NG modes 
arise in /i-SB described by tt/, Hq for /io, and Hq for 
PLh, which we call the phason, the oscillaton, and the 
precesson, respectively (see Table El. The dynamics of 
these three NG modes are similar to those of the phonon 
in the scalar BEG, the magnon in the AFM, and the 
magnon in the EM, respectively. 

The phason is the type-1 NG mode which arises from 
the generator I and similar to the phonon in the scalar 
BEG which describes density fluctuations. Although 
both the oscillaton and the precesson arise from the real 
and imaginary parts, and E^, of the raising operators 
E±a., they are different in their dispersion relation and 
dynamics. These differences arise from the expectation 
value of the Noether charge {Qh^) associated with the 
generator For /xq, both the expectation value (Qh^) 
and the first-order term in dt in the quadratic part of an 
effective Lagrangian in Eq. (I63p vanish. The quadratic 
part of the effective Lagrangian of the oscillaton associ¬ 
ated with the generator E^'^ can be given from Eq. dSB) 
as 

= bo^iVn^f + b^idtTT^f + 5«(V^i)2 + bc^idtirif. 

(77) 

From the effective Lagrangian, two fields and as¬ 
sociated with the generators E^ and E^ are decoupled, 
producing two independent harmonic oscillations of the 
generalized magnetization S'q,: 

( sin[A sin(fc • x)] cos (j) \ 
sin[A sin(fc • a;)] sin(() 
cos[Asin(fe • x)] J 

for ()) = 0 or (78) 

where A (<C 1) represents the amplitude of the oscilla- 
tons, k is the wave number of the oscillatons, and x is 
the coordinate in space. These modes are reminiscent of 
magnons in the AFM, where two magnons representing 
harmonic oscillations of the magnetization appear. On 
the other hand, for fin, neither the expectation value 
{Qhc) nor the first-order term in dt in the quadratic 
part of the effective Lagrangian in Eq. (l63l) vanishes. The 
quadratic part of effective Lagrangian of the precesson as¬ 
sociated with the generator E^’^ is given from Eq. (EOl) 
as 

= p„(n„5tn^ - n:,5tn„) + 6„|vn„p 

= 2ip„(7ri9t7r^ - E^dtiri) 

+6„(V^^)2 + 6„(V7r^)2. (79) 

From the effective Lagrangian, we find that two fields tt^ 
and tTq associated with the generators E^ and E^ form 
a canonical conjugate pair, producing a precession mode 


of the generalized magnetization Sa'. 


Scix) 


( sin A' cos{k ■ x) 
sin A' sin(fe • x) 
cos A' 


(80) 


where A'(<C 1) represents a precession angle. This mode 
is reminiscent of a magnon in the FM, where one magnon 
representing the precession of the magnetization appears. 


E. Differences among three types of NG modes in 
decay processes 

Although both the phason and the oscillaton have lin¬ 
ear dispersions, they play distinct roles in decay pro¬ 
cesses. To see this, let us consider /r-SB with ODLRO 
and fio . Similarly to the proof of Theorem [2 under the 
unitary transformation generated by i?s G l)(s S R’’), 
tt/jIIq, and Ila are transformed into tt/, and 

respectively. Up to the third order in tt/, tt^, 
and TT^, the interaction part of the Lagrangian between 
NG modes, /lint, which is invariant under the transfor¬ 
mations generated by S f) for any s G R’’, should be 
a linear combination of the following four terms: 

( 7 r 7 )^, 7 r 7 n);na for a G 

n^n;n„, n/ 3 n..yn^ for a = /3 -b 7 . ( 8 I) 

Hence we obtain 

/lint = 0(777)^ +777 ^ CaH^Ha 

+ ^ Re(c„,/ 3 ,^n^n;nc) ,(82) 

a,/3,-)'Gfl{p°,a=/3+7 

where c and Ca are real constants, Ca^/s^-y are complex 
numbers, and Re(a:) denotes the real part of x. Thus, 
one phason decays into two phasons tti or two oscillatons 
with the same root vector, whereas one oscillaton Hq 
decays into one phason tt[ and one oscillaton Hq, or two 
oscillatons, lip and H..^ with a = /3 + 7 . 

Similarly, for /x-SB with ODLRO and /X 77 , the interac¬ 
tion Lagrangian between NG modes, /lint, can be written 
up to the third order in 777 , 77 ^ and 77 ^ as 


V. HOMOTOPY GROUPS OF TOPOLOGICAL 
EXCITATIONS FOR p-SB 


Let us now calculate the homotopy groups for /i-SB 
states to find their topological excitations. An element 
of the first homotopy group tti{G/H) characterizes the 
topological charge of a vortex and that of the second 
homotopy group 772 (G/i/) characterizes the topological 
charge of a point defect and that of a skyrmion . Usu¬ 
ally, the homotopy groups are calculated separately for 
individual cases. We here show that not only the homo¬ 
topy groups but also the textures of topological excita¬ 
tions can be calculated systematically for ^-SB. 

We first briefly review the theory of an integral lattice 
and a co-root lattice [4l|, as it is needed for the calcu¬ 
lation of the second homotopy group. We define integral 
lattices Lq and L 77 for a compact Lie group G and its 
subgroup H and a lattice L_r( 5') for a subset S of R+ as 
follows: 


Lq := {/ G R 
L 77 := {/ G R 

Luis) := Spang 


exp(277zi/t) = e,HtG fl}, 
exp(2777i/t) = e,HtG (]}, 

a G S 


(a, a) 


(84) 

(85) 

( 86 ) 


where e is the identity element of G, Ht for t G W is 
defined by 


= (87) 

6=1 


g and () are the Lie algebras of G and H, respectively, 
and SpangA denotes a vector space spanned by elements 
of X with integer coefficients: 

{ m 

y^nkXk 

fe=i 

( 88 ) 



The vector 2cx/ (a, a) and Lfi{R^) are called the co-root 
vector (the inverse root vector) of a root vector a and the 
co-root lattice (the inverse root lattice) of G, respectively. 

Under the addition of r-dimensional vectors, Lq, Lh, 
and L{i{S) {S C R+) form Abelian groups. Let us ex¬ 
amine this point by discussing an example of Lq- Let t 
and s be elements of Lq- The sum t + s satisfies 


Gint = 0(777)^ + 7 ^/ X! CaBj^Ha 
cuGRh 

+ ^ Re (c„,/ 3 ,-^n;n;n„), (83) 

a,(3,‘^GRH ,ck.—(3-{-'y 

where c and Ca are real constants, and Ca,p,-y are com¬ 
plex numbers. Thus, one precesson Hq, decays into two 
phasons ttj and a precesson Hq, or two precessons, H^ 
and H-y with a = /3 + 7 . 


exp(27r7i/t.|_s) = exp(2TriHt) exp{2TTiHs) 

= ee = e, (89) 

and hence we have t + s G Lq- The identity element 
of Lg is the zero vector 0 and the inverse element of 
t G Lg is —t. Since L}i{S) is an Abelian group generated 
by a subset S of R+, L}i{S) is an Abelian subgroup of 
the co-root lattice Lii{R+). Therefore, the coset space 
L[i{R^)/Lii{S) becomes an Abelian group. Writing the 
element of the coset space Lr(Rj^)/Lr[S) as [i] (t G 





9 


TABLE II. Three types of NG modes in /r-SB. The second column shows the classification of /r-SB. The third, fourth, fifth, and 
sixth columns show the generators, fields, dynamics, and dispersion relations of the NG modes, respectively, where tt/ is the 
field of the NG mode generated by the broken generators /, IIq; is defined as Ha := tt^ + iTTa with and being the fields 
of NG modes generated by the broken generators and respectively, and Sa is the generalized magnetization defined 

in Eq. (I42|l . The last column shows a pair of the NG modes into which each of three types of NG modes, 7r/,nQ;(a € i?o) and 
nci;(a G Rh) decay through the interaction Lagrangians in Eqs. (I82II and (I83II . 


name 

classification 

generator 

field 

dynamics 

dispersion 

decay processes 

phason 

ODLRO with fjLH 
ODLRO with /xo 

7 

7 VI 

density 

fluctuation 

linear 

two phasons tt/ 

two oscillatons na(Q: G Rq) 

two precessons na(Q; G Rh) 

oscillaton 

DLRO with /Xo 
ODLRO with /Xo 

Ea and 
(CK ^ Rq) 

(q: ^ Rq) 

oscillation of 

Sa (Eq. (1750 

linear 

one phason tt/ and one oscillaton Ila 
two oscillatons n,g and 11-,, (a = f3 + ~f) 

precesson 

DLRO with fj,H 
ODLRO with fiH 

Ea and E^ 
(a G Rh) 

(a G Rh) 

precession of 
Sc (Eq. ®) 

quadratic 

one phason tt/ and one precesson Ha 
two precessons Ylg and 11-,. (a = /3 -|- 7 ) 


Lr{Rj^), the sum in the coset space Lii{R+)/L r{S) is 
defined as follows: 


[f] + [s] [t + s]. (90) 

It is known that any co-root vector 2Q:/(a,a) is an el¬ 
ement of Lq [m, The co-root vector 2a/(a,a) 

for a S R^\Rh is an element of Lh because G [} 
from TheoremlU where R+\R is the set of elements that 
belong to R+ but not to R. Therefore, 
is a common subgroup of Lji{R^) and Lh, and hence 
the coset space [Lh H Lh{Rj^)\/Lh{Rj\Rh) forms an 
Abelian group. 

Then, the following theorem holds. 


Theorem 3 (1) Provided that iL is a connected sub¬ 
group of G, the first homotopy group for /r-SB is given 
as follows: 


MG/H) 


0 for DLRO; 

Z; for ODLRO and ^ih', (91) 

Z for ODLRO and /.to, 


follows: 


^ 2 (G/iL) = 


Lr(R+) 

Lr(R+\Rh) 

Lr{R+) 

LR(R+\Ri^) 

LHf^LR{R+) 

Lr{R+\Rh) 

LRiR+) 

Lr{R+\RI^^) 


for DLRO and fin', 
for DLRO and /xq; 
for ODLRO and fin', 
for ODLRO and /xq. 


(93) 


Let 6 G [0, tt] and </> G [0, 27 r] be the three-dimensional 
polar coordinates surrounding the point defect, 0{9, (j)) 
be the value of the order parameter at {9, </>). The texture 
of the point defect with topological charge [ 2 a/(a, a)] is 
given by 


0(9, (j)) = exp 




(a.,a)_ 


o exp 


i9 


Ei 


(a,a)_ 


O On 


(94) 


where Oo and g o Oq (g G G) is the value of the order 
parameter at 0 = ^ = 0 and that of the order parameter 
obtained by the symmetry transformation g from Oq, 
respectively. 


where Z; = {g^\i = 0 , 1 , 2 , • • • , Z — 1 } is the cyclic group 
of order I (g is the generator of Z/), Z is a positive inte¬ 
ger which is uniquely determined from fin, and Z/=i is 
a trivial group consisting of the identity element alone. 
Let 9 G [0,27r] and 0(9) be the azimuth angle around 
the vortex and the value of the order parameter at the 
angle 9. The vortex with topological charge g represents 
a vortex around which the phase of the order parameter 
rotates by 27r: 

0(9) = e*^Oo, (92) 


The proof of this theorem is given in Appendix □ 

Physically, nontrivial elements Z; and Z in Eq. m 
represent quantum vortices. We note that a quantum 
vortex does not necessarily have an integer charge in a 
system with internal degrees of freedom [43l |. Let us ex¬ 
amine this point as yet another example of /x-SB. Con¬ 
sider the three-dimensional representation of g = su(N = 
2) and a phase with ODLRO and fin- This example will 
appear as a special case of = 2 of a Z7(iV)-symmetric 
system in Sec. I VII In this representation, the order pa¬ 
rameter (As) is a 2 X 2 symmetric complex matrix which 
transforms under the symmetry transformation of U(2) 
as follows [13 : 


where Oq is the value of the order parameter at 0 = 0 . 
(2) The second homotopy group for /x-SB is isomorphic 
to the coset space constructed from the co-root lattice as 


(As) C/(As)‘Z7 for U G U(2). (95) 

In /i-SB with ODLRO and ^h, the expectation value of 
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(As) is given by 

(As) = (Ao) = ( J S ) • (96) 

Since this representation is the three-dimensional rep¬ 
resentation of U{2), which is equivalent to the three- 
dimensional representation of U{1) x SO{3), the order 
parameter (As) is related to the three-dimensional com¬ 
plex vector {(f>) = (((/)i)' ((t>f)), (d)-i)), which is the order 
parameter of a spin-1 BEC [111, [13; as 

= ( (« 

The first homotopy group tti (G/H) of this phase is given 
by 

Tri{G/H) = Z 2 = {e,g}, (98) 


where e and g are the identity element and the gener¬ 
ator of Z 2 with g^ = e. From Eq. (EH), a vortex with 
topological charge g is given by 


(Af)(0) = 




0 

0 0 


(99) 


where 6 * G [0,27r] represents the azimuth angle around the 
vortex. We can show = e as follows. Let CFx,<Jy,az 
be the Pauli matrices. From Eq. (EH, the texture of the 
vortex with topological charge g^ is given by 


(Af m = 

U{9) = e 


^2ie Q 

0 0 


= Ui0){AoyU{9), ( 100 ) 


iBcTz 


( 101 ) 


By the continuous deformation defined by 


(Asy)(0) =t/(0,t)(Ao)‘t/(0,t), 

U{e,t) =e-”‘^e*®^e” 




( 102 ) 

for 0 < t < 1 , 

(103) 


(As^t) is transformed from (Ag^t^o) = (Af ) to a uni¬ 
form order (Asy=i) = (Aq), which implies = e. 

From Eqs. (15^ and EH, 7 r 2 (G/iL) is generated by 
a set of representative elements {[ 20 ./{a^<y)\\a. G i?+} 
of the coset space. This shows that any point defect 
in /r-SB can be represented as a composite of the point 
defects with topological charge \2ot/{a.,a.)\ {a G R+)- In 
Eq. dHH, the coset spaces of Lr{Rj^) by their subgroups 
Lr{S) [S — R-\.\Rh T R+\Rq T R+\R^^) are considered 
instead of the numerators Lr{R+) or L//nLi{(i?+). This 
is because the element of the denominators Lfi{S) of the 
coset does not give nontrivial topological excitations. Let 
us clarify this point for the case of iiR- Since and 
Ha for a G R+\Rh are unbroken generators for fin, 


the successive action of exp 


Ha, 


{CK.,Ct) 


and exp 


ie 


{ck.,cx) 


leaves the order parameter invariant: 


exp 


Ha 


(a, a) 


o exp 


ie- 


E. 


•I ^ 


O On = On 


(a,a)_ 

forV(0, (^) G [0, tt] X [0, 27r]. (104) 


Therefore, Eq. (IMl) does not give a nontrivial point defect 
but a uniform order for a G R+\Rh for ^ih- 

The point defect in Eq. (1941) with topological charge 
2Q:/(a, a) is similar to the point defect in the EM and the 
AFM in that the former is obtained by replacing su(2) 
spin S by the generalized magnetization Sa- In fact, 
let M be the magnetization of the EM or that of the 
sublattice in the AFM. A point defect can be described 
as a hedgehog configuration of M: 

M{9,(j)) = (sin0cos(/), sin^sin^, COS0), (105) 

where 9 and cj) are 3-dimensional polar coordinates. This 
hedgehog configuration is obtained by the successive ro¬ 
tation of the spin around the y-axis by angle 9 followed 
by the rotation around the z-axis by angle (f>: 

M{9, (j)) = exp {i(j)Sz) o exp (i9Sy) o Mq, (106) 

where Mq = (0,0,1) is the magnetization of the EM or 
that of the sublattice in the AFM parallel to the z-axis. 
Comparing with Eq. (PI). we can see that the point defect 
in Eq. (1941) is the generalization of the point defect in the 
EM and the AFM obtained by replacing S hy Sa- 


VI. APPLICATION OF g-SB TO 
[/(A)-SYMMETRIC systems 

In this section, we apply g-SB to {7(V)-symmetric sys¬ 
tems [Il[Ii,|4llil. Here, we consider up to the third 
lowest-dimensional representation of the su(V)-Lie alge¬ 
bra since symmetry broken phases are not necessarily 
characterized by g-SB in higher-dimensional representa¬ 
tions. Up the third lowest-dimensional representation, 
the irreducible representations of the su(V)-Lie algebra 
are given by the following three representations, the N- 
dimensional representation, the N{N — 1)/2-dimensional 
representation, and the N{N -|- l)/2-dimensional repre¬ 
sentation. To make this paper self-contained, we briefly 
review the symmetry transformation, the highest weight, 
and the set of wei ght vectors of these three represen¬ 
tations. See Refs. |20l . for detail on the su(fV)-Lie 
algebra and its representations. 

Let be a set of real (iV—l)-dimensional vectors 

that satisfy 

N 

'£i^j=0, (107) 

J=1 

=6ij - (108) 

Defining cti j by ctij := Ui — Uj, the set of positive root 
vectors i?+ of the su(V)-Lie algebra is given by 

R+ = {aijlotij := Ui - <i < j < N}. (109) 

The weight vector Vi is normalized in Eq. (11081) so that 
the magnitude of the root vectors is \/2. For example. 
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for TV = 2 and 3 are given by 

N = 2:., = ±,., = -±, ( 110 ) 

.3=(0.^). (Ill) 

For N = 2 and iV = 3, the schematic illustrations are 
presented in Fig. [T](a) and (b), respectively. 

The lowest-dimensional representation is the N- 
dimensional representation. The element of this repre¬ 
sentation is an A^-dimensional complex vector v. The 
symmetry transformation of this representation acts on 
V as an action of the matrix from the left: 

V 1 -^ Uv ioT U G U{N). (112) 

For the ODLRO, v is the order parameter of a BEC in 
degenerate A^-component bosons. The highest weight fin 
and the set of weight vectors W[fiH] are given by 


/x// = ivi, (113) 

W[u,] = {u,\j = l,2,--- ,N}. (114) 

The weight vector state |z^j) is a unit vector whose com¬ 
ponents vanish except for the j-th component: 

|«^,)=*(0,...,0,lo,...,0). (115) 


The antisymmetric nature of arises from the anti¬ 
commutation relation of the fermions. Under the U{N)- 
symmetric transformation, ipi and transform as 

!-)■ (119) 

( a ^)^ ^ U,k{A^Pi)CU)i, = (uAju)^ 

for U e UiN), (120) 

which coincides with the transformation ()116l) . The high¬ 
est weight fiH and the set of weight vectors W[fiH] are 
given by 


flH = 1^1 + 1^2, ( 121 ) 

-I- v>2] = {i^i + < i < j < N}. (122) 


The weight vector state -|- i^j) is the skew symmetric 
matrix whose elements vanish except for the (i, j)-th and 
(j, i)-th elements: 


■^) = 





(123) 

(124) 


The third lowest-dimensional representation is the 
N{N + l)/2-dimensional representation. The element 
of this representation is an x TV complex symmetric 
matrix Ag. In fact, the dimension of the set of TV x TV 
complex symmetric matrices is TV(TV -|- l)/2. Similarly to 
Eq. (11161) , the symmetry transformation of this represen¬ 
tation acts on Ag as 


Ag UAs*U for U e U{N). (125) 


The second lowest-dimensional representation is the 
TV (TV — 1)/2-dimensional representation. The element of 
this representation is an TV x TV complex skew symmetric 
matrix A^. In fact, the dimension of the set of TV x TV 
complex skew symmetric matrices is TV(TV — l)/2. The 
symmetry transformation of this representation acts on 
Aa as an action of the matrix and that of its transpose 
from the left and the right, respectively: 

Aa UAjU for U e UiN). (II6) 

For the ODLRO, this order parameter corresponds to 
that of an s-wave superfuild phase in degenerate TV- 
component fermions in a non-relativistic system [23 - 1^ . 
In fact, let be the fields of the degenerate TV- 

component fermions. Since the TV-components are de¬ 
generate and the total number of fermion is conserved in 
a non-relativistic system, this system is invariant under 
the U (TV)-symmetry transformation: 

t/i '->■ Uijipj for U G UiN). (117) 

Here, repeated indices are assumed to be summed over 
i = 1, - ■ ■ , TV. The order parameter of the phase is given 
by the following TV x TV matrix: 

Aa = ( 118 ) 


Although the transformation in Eq. (I125|) coincides with 
Eq. (I116p . Aa has the odd parity, ‘Aa = — Aa, and 
Ag the even parity, ‘Ag = Ag. An TV x TV skew sym¬ 
metric (symmetric) matrix A transforms into a skew 
symmetric (symmetric) matrix under the transformation 
A I— UA*U ([/ G UiN)). For the ODLRO, the order 
parameter Ag is related to that of a p-wave superfuild 
phase in degenerate TV-component fermions [l^. In fact, 
let {V'iliLi be the fields of the degenerate TV-component 
fermions. The order parameter of the p-wave superfluid 
phase is given as follows (47l| : 

{'4^i,h'4^j, — k) — ^ ( ka^Aai ij. (126) 

OL—x,y,z 

where three TV x TV matrices A^ := (a = 

x,y,z) are symmetric matrices (47l| . The symmetric na¬ 
ture of Aa arises from the anticommutation relation of 
the fermions and the odd parity of the orbital part of 
the p-wave pairing. Under the symmetry transformation 
that mixes the degenerate TV-components, ipi i—>■ Uijipj, 
the N X N symmetric matrix Ag transforms according to 
Eq. (11251) . The highest weight fin and the set of weight 
vectors Wlfin] are given by 

Uh = 21^1, (127) 

W[2iyi] = {u, + Uj\l<i<j < N}. (128) 
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The weight vector state \ui + 1^^) is the symmetric matrix 
whose elements vanish except for the (i, j)-th and (j, *)-th 
elements: 


and hence is a connected group. For the case of /xq, 
from Eq. (15^ the expectation value of the Casimir in¬ 
variant within the unit cell u is given by 


A. 



(129) 


i^u 

^(77) 

kl 

(130) 

G| = 


(137) 


Classification of fi-SJi phases 


From Eq. (11071) , the right-hand side of Eq. (I137p vanishes 
when the unit cell consists of iV-different weight vectors, 
Thus, the ground state is given by 


We first classify the /r-SB phases for g = su{N) that 
appear up to the third lowest-dimensional representation. 

1. Lowest-dimensional representation 

For /r-SB with ODLRO and fin, the expectation value 
of the order parameter coincides with the highest weight 
of this representation in Eq. (11141) : 

(0) = |ivi) = (l,O,--- ,0). (131) 

The remaining symmetry H of the state is given by 

H :={U €U{N)\U\ui) = livi)} 

= {U€ U{N)\U,, = [/,! = (5,1 (j = 1,2, •. ■ , N)} 
-U{N-1), (132) 

where ~ represents the group isomorphism. Thus, 
is a connected group. In this representation, the 
Casimir invariant C'|((</>)) coincides with \{4>)\‘^ and 
hence Cf ((</>)) ^ 0 so long as the order parameter has a 
nonzero expectation value. Therefore the pair of ODLRO 
and /xq is absent in this representation. 

We next consider the case of DLRO. There are A^-fold 
degenerate states in each site which are referred to as 
color or flavor. For /r-SB with DLRO and /jlh, the mean- 
field ground state is given by 


|G5) = (g)(g)lM.)^, (138) 

u^lA iGu 

{l^i}ieu = {I'j} j^i- (139) 

This is the mean-field ground state of the so-called N- 
color density wave (iV-CDW) phase [48l - l5lj| . This state is 
a generalization of the S't/(2)-antiferromagnet to a gen¬ 
eral SU{N)-spm system. In the former, two states, a 
spin-up state and a spin-down state, constitute a unit 
cell, while in the latter N states do. The unitary op¬ 
erator U associated with the symmetry transformation 
U G U(N) leaves [G^) unchanged up to a global phase 
factor if and only if U leaves all of the iV-different weight 
vectors |iZi) within each unit cell unchanged up to a global 
phase. In other words, the condition 

C7|iz,) = e*'^|iz,) for Bcj) G R, (140) 

must be satisfied for all of TV-weight vectors 
This condition is satisfied if and only if C/ is a diagonal 
matrix. Thus, the remaining symmetry H of the state is 
given by 

iL :={C/Gt/(iV)|C/ = diag(e*■^^e*‘^^••• ,e*'^"), 

Uaf, (141) 

and hence 77 is a connected group. 


|G5) =(g)|i^i)*. (133) 

i^L 

This is the mean-field of an 5'7/(A^)-ferromagnet (43| . 
The unitary operator U associated with the symmetry 
transformation U G U{N) leaves IGS”) unchanged up to 
a global phase factor if and only if U leaves the weight 
vector \ui) on each sites unchanged up to a global phase: 

[/|jzi)=e*‘^|izi)for3(;5GM. (134) 

This condition is equivalent to the condition that the first 
row and the first column of U vanish except for a diagonal 
element: 

Uij = Uji = 0 for j = 2,3, • • • , TV. (135) 

The remaining symmetry 77 of the state is given by 

77 := {7/ G C/(iV)|[/i, = Uji = 0 (j = 2, 3, • • • , TV)} 

~ 77(1) X 7/(A^-1), (136) 


2. Second lowest-dimensional representation 

For /i-SB with ODLRO and fiH, the expectation value 
of the order parameter is given from Eqs. ([8]) and (1531) by 

Aa=|«^i+i^2) = Aii’2)_ (^42) 

The remaining symmetry 77 of the state is determined 
by a straightforward calculation as 

77 := {77 G 77(iV)|77Aii’2)*77 = A^i’^)} 

= (77 G 7/(iV)|77ii7722 - U 12 U 21 = 1, 

7/y =77,,=0(z = l,2,j = 3,4,..-,iV)} 

~ 57/(2) X 77(iV-2), (143) 

and hence 77 is a connected group. On the other hand, /i- 
SB with ODLRO and ^jlq does not necessarily exist for ar¬ 
bitrary N. The same point is discussed in the derivation 
of the mean fields of /r-SB in Sec. [TTl The set W[ui -\- 1 / 2 ] 
does not necessarily contain the zero-weight vector (^ . 
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In fact, Vi + Vj = 0 {1 < i < j < N) only when the fol¬ 
lowing conditions are met: N = 2, Vi = Vi, and Vj — 1 ^ 2 ■ 
For N greater than 2, W\vi + 1 / 2 ] does not contain the 
zero-weight vector. 

We next consider the case of DLRO. There are N{N — 
l)/2-degenerate states labeled by a weight vector fi S 
W[vi + V 2 ] in each site. For /r-SB with DLRO and fin, 
the mean-field ground state is given from Eqs. m and 
(HSl) by 


\GS)=<^\v,+V2)^. (144) 

i&L 

Similarly to the DLRO and fin in the lowest-dimensional 
representation, the unitary operator U associated with 
the symmetry transformation U G U{N) leaves \GS) un¬ 
changed up to a global phase factor if and only if U leaves 
the weight vector \vi+V 2 ) on each site unchanged up to 
a global phase: 

U\vi + V 2 ) = e*‘^|i>'i -f V 2 ) for G R 

=C/,i=0 (^ = l,2,j = 3,4,••• ,Af). (145) 

Thus, the remaining symmetry H of the state is given by 

H := {U G U{N)\ 

t/,=0 (z = l,2,j = 3,4,-.- ,iV)} 

~ C/(2) X C7(iV - 2), (146) 

and hence H is a connected group. For the case of fiQ, 
the right-hand side of Eq. (11371) vanishes when the set of 
weight vectors {/XijiguOn the unit cell u is 

I = {V 2 j-i + for even TV; 

1 = {vj + U {vi + vn} for odd N. 

(147) 

For even N, TV/2-sites within each unit cell are sufficient 
because we have from Eq. (11071) 

N/2 N 

+ V2j) = ^ IZj = 0. (148) 

On the other hand, (TV — l)/2-sites within the unit cell 
are not sufficient because the sum of the weight vectors 
within each set i® 

(JV-l)/2 AT-i 

{V2j-1 + U2j) = ^ 

1=1 1=1 

We have to consider the unit cell with iV-sites. In fact, 
the set {vj + U {vi + izjv} satisfies 

N-l N 

Y ('^T + '^T+i) + (*^1 +vn) = 2Y^j=^ (150) 
1=1 1=1 


Thus, the ground state is given by 

|G5) = (g)(g)|/.t0„ 

u^lA i^u 

[ {Milieu = {V2j-l + V2j}^ll 
\ (Milieu = 1*^1 + *^l + l|jLl^ O {ui -|- IZtvI 


(151) 

for even TV; 
for odd TV. 

(152) 


The remaining symmetry H of the state can be calculated 
in a manner similar to the derivation of Eq. (11411) . A 
unitary matrix U is included in H if and only if U\fii) = 
e*'^|Mi) G 1^) for all of the weight vectors within each 
unit cell. For even TV, U\v 2 j-i -I- i>' 2 j) = e^^\v 2 j-i -I- 
V2j) {3(j) G M) is satisfied when both the {2j — 1) and the 
(2j)-th columns and the (2j —1) and the (2j)-th rows van¬ 
ish except for the (2j - 1, 2j - 1), {2j - 1, 2j), (2j, 2j - 1), 
and (2j, 2j)-elements. Therefore, U is block-diagonalized 
into a direct product of 2 x 2 matrices and we obtain 

U, ([/, G U{2)) I , (153) 

~ [^7(2)]^/^ (154) 

For odd TV and the set of weight vectors {uj + Vj+i}Y^^U 
(izi -|- izatI, a unitary matrix U is included in H if and 
only if [/ is a diagonal matrix. Therefore, we obtain 

iL :={Dg D(TV)|D = diag(e*'^Se*'^=,--- ,e*‘^'') (/>. G K| 
Uilf, (155) 



For both even and odd TV, iV is a connected group. 


3. Third lowest-dimensional representation 

For fi-SB with ODLRO and fin, the expectation value 
of the order parameter is given by 

A« = |2izi) = All’ll. (156) 

The remaining symmetry H of the state is given by 
H — {Ug C/(TV)|C/a 11’11*C/ = Ali’i)} 

= {Ug U{N)\U,, = Uji = 0 (j = 1,2, • • • , TV)} 
~C/(TV-1), (157) 

and hence H is a connected group. Similar to the case of 
the second lowest-dimensional representation, /r-SB with 
ODLRO and /xg does not necessarily exist for arbitrary 
TV. The set VF[2i>'i] contains the weight vector Vi Vj = 
0(1<*<J< TV) only when the following conditions 
are met: N = 2, Vi = Vi and Vj = U 2 [2^. For TV 
greater than 2, VF[2izi] does not contain the zero-weight 
vector. For the ODLRO and /xg with TV = 2, 77 is not 
a connected group. In fact, the order parameter of this 
phase and H are given by 

As = (Jx, (158) 

7J={e*‘/2|tGR| X {/2,e*^CTy} 

= 7/(1) X Z 2 , 


( 159 ) 
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where TV x TV' is a semidirect product whose product is 
given by 

{n,h) X [n',h') = {nhn'h~^,hh') 

for Vn, n' G N, V/i, h' gN'. (160) 

We next consider the case of DLRO. There are N{N + 
l)/2 states labeled by a weight vector /x G Wpixi] in 
each site. For /x-SB with DLRO and ^jlh, the mean-field 
ground state is given by 

|G5) = (g)|2ixi),. (161) 

i&L 

The remaining symmetry H of the state is calculated in 
a manner similar to Eq. (11361) as 

H:={Ug UiN)\UAi^’^^*U = for Bcj) G K} 

= {Ug UiN)\Uij = Uji = 0 (j = 2,3, • • • , N)} 

- U{l)xU{N-l), (162) 

and hence iJ is a connected group. For the case of /xq, 
the ground state and the remaining symmetry H of the 
state is obtained in a manner similar to the case of the 
DLRO and /xq for the lowest-dimensional representation, 
the right-hand side of Eq. (11371) vanishes when the set of 
weight vectors {/Xi}ig„ within the unit cell u are 

{/x*Kg. = (163) 

Thus, the ground state is given by 

|G5) =(g)(g)|/x.)^, (164) 

u£U i^u 

{/x,},e. = (165) 

The remaining symmetry H of the state can be calculated 

in a manner similar to the derivation of Eq. (11411) . A 

unitary matrix U is included in H if and only if [/ is a 
diagonal matrix. Thus, the remaining symmetry H of 
the state is given by 

iL :={t/e[/(^^)|^7 = diag(e*'^^e*'^^••• ,e*'^'^), <()* G K} 
~ t/(l)^, (166) 

and hence iL is a connected group. 


B. Numbers of NG modes 

We next calculate the number of NG modes in /x-SB 
phases classified in the previous subsection. From the 
quadratic part of effective Lagrangians in Eq. (15^ . the 
numbers of type-1 and type-2 NG modes, ni and 7X2, are 
given by 

for DLRO and fin; 
for DLRO and /xq; 
for ODLRO and /x/x; 

[ {2\R^^\ -f 1,0) for ODLRO and /xq, 

(167) 

where |A| denotes the number of elements in the set X. 


(ni,n2) = 


(OARhI) 
(2|i??|,0) 
(1> \Rh\) 


1. Lowest-dimensional representation 

Combining the set VF[/xjx] of the weight vectors in 
Eq. (11141) and the set of the weight vectors in 

the unit cell u in Eq. (I139|) . the sets Rh and R^ are 
given from Eqs. ((171) and ((15|) by 

i?xx = {q;ij |2 < j < iV}, (168) 

R^ = {a,^j\l<i<j<N} = R+. (169) 

Substituting these equations into Eq. (11671) . we obtain 

( {0,N — 1) for DLRO and /xjx; 
( 711 , 7 x 2 ) = < (iV(iV —1),0) for DLRO and/X q; 

y ll,N — 1) for ODLRO and /xjx. 

(170) 


2. Second lowest-dimensional representation 


Combining the set VF[/xjx] of the weight vectors in 
Eq. (11221) and the set {/Xi}ig„ of the weight vectors in 
the unit cell u in Eq. (I152|) . the sets Rh and i?g are 
given from Eqs. ((471) and ([48]) by 

Rh = K = 1,2, 3 < j < N}, (171) 

i?o = R-e\{oi 2 j-i, 2 j\j = 1,2,--- ,fV/2} for even N, 

(172) 

i?Q = {clij |1 < i < j < N} = i?+ for odd N. (173) 


For N = 2, there exists the ^-SB phase with ODLRO and 
/Xq. From Eq. (|4^ . the set i?Q ° in this phase is empty: 

= 0, (174) 


where 0 denotes the empty set. Substituting the above 
three equations into Eq. (11671) . we obtain 


(^ 1 , 7 x 2 ) 

r (0,2(fV-2)) 
{N{N-2),0) 
= { (fV(fV-l),0) 
(l,2(fV-2)) 

I (1,0) 


for DLRO and /x/x; 
for DLRO, /xo,and even N; 
for DLRO, /xo,and odd N; 
for ODLRO and /xxr; 
for ODLRO, /xg, and N = 2. 

(175) 


3. Third lowest-dimensional representation 

Combining the set VF[/xjx] of the weight vectors in 
Eq. (I128P and the set {/Xi}ig„ of the weight vectors in 
the unit cell u in Eq. (I165() . the sets Rr and i?g are 
given from Eqs. ((471) and ([48]) by 

Rh = {ai,x |2 < j < N}, 

Rq = {ajj |l < X < j < A^} = R+- 


(176) 

(177) 
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For N = 2, there exists a /i-SB phase with ODLRO and 
fiQ- The set in this phase is given from the definition 
of in Eq. (|4^ as 

= {ai.2}. (178) 

Substituting the above equations into Eq. (11671) . we ob¬ 
tain 


{ni,n2) 

I (iV(iV-l),0) 
] ( 1 ,^- 1 ) 

I (3,0) 


for DLRO and ^ih\ 

for DLRO and /.tg; 

for ODLRO and 

for ODLRO, /xo,and N = 2. 

(179) 


For the ODLRO and it is easier to calculate 
tti{G/H) and 712 (G/i?) directly rather than using The¬ 
orem [3] because the order parameter manifold is isomor¬ 
phic to a higher-dimensional sphere: 

G/H = U{N)/U{N -l) = (185) 

Since TV > 2, we obtain the following two equations from 
the standard results of homotopy groups [4^ : 

^i(G/i7) = 7ri(52^-i) =0, (186) 

7r2(G/i7) = 712 ( 52 ^-^) = 0. (187) 


2. Second lowest-dimensional representation 
For the DLRO, from Eq. diu, we obtain 


C. Homotopy groups of topological excitations 

Finally, we calculate the first and second homotopy 
groups for /r-SB phases. Since H is a connected group in 
all the cases except for the ODLRO and fiQ with N = 2 
in the third lowest-dimensional representation, we can 
apply Theorem [3] except for this case. 


7ri(G/i7) = 0. (188) 


For the DLRO and fin, from Theorem [3] and Eq. (I173L 
we obtain 


TT2{G/H) = 


Lr{R+) 

T_R(i?+\i?ff) 


= Spanz{Q;2,3} - 


(189) 


For the DLRO, /xg and even TV, from Theorem |3] and 
Eq. (j 15211 . we obtain 


1. Lowest-dimensional representation 

For the DLRO, from Eq. (IMl) . we obtain 
^i(G/i7) =0. 

For the DLRO and fin, we obtain from i?+ in Eq. (11091) 
Lr{R+) = Spanz ■ 


Lr{R+\Rh) = Spanz {a,,, |2 < i < TV} , (182) 




( JV/2 

] 


L«(i?+\i?OD) = 

} ^ ]'rnja2j-i,2j 
\ ? —1 

ruj G Z > , 

(180) 



) 

(190) 

(fTMl) 

T^2{G/H) = 

Lr{R+) 


Lr{R+\RO^) 



f 2a 

cx. € R-\- ^ ^ 

" Ar/2-l 


1 (a, a) 

yy mja2j,2j-ei 

ruj G Z 

VI 

i<N}. (181) 

1 ^=1 

J 




(191) 


In deriving the first equality of Eq. (|181L we use 
Eq. (11091) and = 2 for any i,j. Therefore 

we obtain from Theorem [3] 


7T2{G/H) 


Lr{R+) 
Lr{R+\Rh) 
Spanz{Q:i,2} — 2- 


(183) 


For the DLRO and /xg, since Rg = R+, we obtain from 
Theorem [3] 


7r2{G/H) = Lr{R+) 


N 


— TnjkOLj,k 

T,fc=i 


mjk G Z 


N 


= \ 

[t=i 


N 


m 


'3 G Z, rrij — 0 
T=i 


( 184 ) 


For the DLRO, /xg and odd TV, we obtain the same 
result Tr 2 {G/H) = as in the case of the lowest- 

dimensional representation since R^ is the same in 
both cases of the lowest-dimensional and second lowest- 
dimensional representations. 

For the ODLRO and fiR, we can prove -KiiG/H) — 0 
as follows. To show this, it is sufficient to show that the 
following vortex-like texture analogous to Eq. (IMI) can 
be deformed into a uniform order: 

(A„)=e*«Aii’2). (192) 

Here, 9 S [0,27r] is the azimuth angle around the vortex¬ 

like object. Let and \y be two matrices defined by 

1 0 0 \ 

0 0 0, (193) 

0 0 -1 y 

0 0 -T\ 

0 0 0, (194) 

T 0 0 y 
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and define N x N matrices and Ay by a direct product 
of A^ and Ay with the identity matrix Ins of size (iV—3), 
respectively: 


Az — Az © Ins, 

\ = ® Ins- 

By using Az, Eq. (11921) can be written as 
(A„) = C/(d)A(i' 2 )it/( 0 ), 

Consider the continuous deformation defined by 


(Aa.t)(0) = [/(0,<)(A(i'2))‘l7(0,t), 
{7(0, t) = 

for 0 < t < 1. 

The unitary matrix U{0,t) satisfies 


(195) 

(196) 


(197) 

(198) 

(199) 

( 200 ) 


For the DLRO and fj-H, we obtain the same result 
7r2 (G/iJ) = Z as in the case of the lowest-dimensional 
representation since Rh coincides in both cases of the 
lowest-dimensional and third lowest-dimensional repre¬ 
sentations. For the DLRO and fj-o, we obtain the 
same result 7r2(G/i7) = Z^“^ as in the case of the 
lowest-dimensional representation since coincides in 
both cases of the lowest-dimensional and third lowest- 
dimensional representations. 

For the ODLRO and fin, we can prove 'KiiG/H) = Z 2 
in a manner similar to the discussion in Sec. Let g be 
the generator of Z 2 . The vortex with topological charge 
g'^ is described as 

As(0) = e2*®A(i’i\ (209) 

where 0 G [0, 27r] is the azimuth angle around the vortex. 
We can prove = e in a manner similar to the case 
of TV = 2 in Sec. |V] by replacing Pauli matrices ai{i = 
X, y, z) into 


{7(0, t = 0) = exp(i0Az) = C/(0), (201) 

{7(0, t = 1) = 

= =/Ar, (202) 

where In denotes the identity matrix of size N. In the 
third line, we use the relation 

e-^’^^Aze*’^^ = -A^. (203) 

Therefore, (A^ j) is transformed from Eq. (119211 to a uni- 

form order {A.a^t=i) = Aa ’ , which implies the triviality 
of the vortex-like texture in Eq. (I192|) . From Eq. (I173|) . 
we obtain 


Lh 



t, rrijk SM, 3<j<j<iV. 


J 


(204) 


Lh n Lfl(i?+) 

= Span^ [{cKi.i|3 < i < j < N} \J {ai, 2 }] 

= LRiR+\RH). (205) 

Therefore, we obtain tt 2 {G/H) = 0 from Eq. (IMl) . For 

the ODLRO and fiQ with iV = 2, by using Eq. (ICT)) and 
substituting 7?°^ = 0 into Eq. (1^ . we obtain 

7ri(G/77)=Z, (206) 

TT2iG/H)= 0. (207) 


3. Third lowest-dimensional representation 

For the DLRO, from Eq. (IMl) we obtain 

7ri(G/i7)= 0. (208) 


■—Hi® lNs{i = x,y,z). ( 210 ) 

For the ODLRO and we obtain the same result 
7r2 (G/i7) = 0 as in the case of the lowest-dimensional 
representation since Lh and Rh coincide in both cases 
of the lowest-dimensional and third lowest-dimensional 
representations. For the ODLRO and /xq with N = 2, 
Eq. m is not applicable because H is not a connected 
group. From the correspondence with the spin-1 BEG 
in Eq. dSil), this phase coincides with the polar phase in 
the spin-1 BEG The first and second homotopy 

groups of this phase are given as follows [H, : 


^i(G/77)=Z, (211) 

tt2{GIH)=Z. ( 212 ) 

We list the results obtained in this section in Table uni 
together with the examples of the classified phases. Ex¬ 
amples with g = u(l)©so(3) are included because so(3) is 
isomorphic to su(2). From the fifth column of Table lllll 
we can see that a large class of symmetry broken phases 
are described in terms of /i-SB. 


VII. DISCUSSION ON THE CASE OF 
HIGHER-DIMENSIONAL REPRESENTATION 

So far we have confined our discussions to low¬ 
dimensional representations. We next turn to the case 
of a higher-dimensional representation. In the case of 
ODLRO in a higher-dimensional representation, there 
appears more than one Gasimir invariant in the en¬ 
ergy functional. Due to the competition between these 
Gasimir invariants, the phases that arise from the min¬ 
imization of the energy functional are, in general, de¬ 
scribed by neither {^-SB nor inert states, where an inert 
state is a state in which the order parameter is indepen¬ 
dent of the coupling constants |5^ . In this section. 
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TABLE III. Classification of /.i-symmetry breaking in systems without Lorentz invariance. The Lie algebras of the systems are 
assumed to take the form of g = u(l) © su{N) with N > 2. Here, is the set of weight vectors in an Ai-dimensional 

representation of su(A^), which satisfies Eqs. (I107II and (I108II . The first column shows the irreducible representation of the order 
parameter (</>) for ODLRO and that of the field of particles on each site of the lattice L for DLRO. The second column shows 
whether the system is characterized by fiH or fxo and by ODLRO or DLRO, respectively. The expectation value of the order 
parameter for ODLRO can be written in the form of Eq. m or Eq. while the ground state for DLRO can be written in the 
form of Eq. GSI or Eq. dJ. The row with /tq and ODLRO appears only for N — 2. For the representation of the /irr = i^i 
case, the pair of /iq and ODLRO is absent since this representation is the lowest-dimensional one. The third column (ni,n 2 ) 
shows the numbers of type-1 and type-2 NG modes. The fourth column (7ri,7r2) lists the first and second homotopy groups 
of G/H. FM, AFM, CDW, SF, and VBS stand for ferromagnet, antiferromagnet, color density wave, superfiuid, and valence 
bond solid, respectively. In the row with 2vi, fio, and DLRO, the upper (lower) row corresponds to the case of even (odd) N. 


representation ij,h classification 

{ni,n2) ( 711 , 712 ) 

example 

izi ODLRO, fj,H 

(1,A-1) (0,0) 

SU{N)-FM BEC 

DLRO, fiH 

(0,A-1) (RZ) 

spin-1/2 FM, SU{N)-FM 

DLRO, po 

(A(A- 1),0) 

spin-1/2 AFM, A-CDW [48-^ 

vi + V2 ODLRO, /xif 

(l,2(A-2)) (0,0) 

s-wave SF in 3-component fermion [24-26] 

ODLRO, /r.0 

(TO) (z;o) 

s-wave SF in 2-component fermion [3J 

DLRO, fiH 

(0,2(A-2)) 


DLRO, fjLQ 

(A(A-2),0) (0,ZT-i) 

VBS in SU (4)-spin model [31., .M] 

(A(A- 1),0) (0,2"''--') 


2v\ ODLRO, fj,H 

(1,A-1) (Z2,0) 

FM phase in spin-1 BEC [21., .22] 

ODLRO, po 

(300) (z;z) 

polar phase in spin-1 BEC 121. 221 

DLRO, fiH 

(0,A-1) (RZ) 

spin-1 FM 

DLRO, fio 

(A(A- 1),0) (0,Z"''--') 

spin-1 AFM 


we focus on the case of a higher-dimensional representa¬ 
tion in which two Casimir invariants appear in the en¬ 
ergy functional. In this case, many of the ground states 
are described by inert states despite the competition be¬ 
tween Casimir invariants. Let us examine this point by 
discussing examples of ^in-2 BECs (TB - [T7l| and spin-1 
color superconductors [13 ■ 


are the highest-weight, the zero-weight, and the lowest- 
weight states, respectively. For g = u(l) © so(3), the 
lowest-weight state is obtained by applying 7r-rotation 
around the x-axis to the highest-weight state. The mean- 
field energy functional can be constructed from the norm 
I (0)1 of the order parameter and the Casimir invariants 
of SO{5) and SO{3) as 


A. Spin-2 BEC 

First, we consider the example of spin-2 BECs. As we 
will see below, all of the ground states are inert states. 
The symmetry group of the system is U (1) x SO{3). The 
order parameter of spin-2 BEC is a five-dimensional com¬ 
plex vector 

(0) = *((02), (0l), (0O), (0-l), (0-2)), (213) 


n(0))=-c|(0)P+c'|(0)|4+c;Cr^^)((0)) 

(218) 

where and ((0)) are the Casimir in¬ 

variants of 50(5) and 50(3), respectively [l^. Here 
02 °^^^((0)) is related to the spin-singlet pair amplitude 


and the Cartan generator of 50(3) is the 5^-operator 
defined by 


5, 


/2 0 0 0 0 \ 

0 10 0 0 
0 0 0 0 0 

0 0 0 -1 0 

\0 0 0 0 - 2 / 


(214) as 


The three eigenstates of the 52-operator 


Aoo(0) -^[020-2 — 010-1 + (0o)^] (219) 


cf^\m = m\'^-5\Aoomr- mo) 


( 0 )^^ =‘( 1 , 0 , 0 , 0 , 0 ), 

(0)2 = ‘(O,O,1,O,O), 

(0)l =‘(0,0,0,0,1), 


(215) 

(216) 
(217) 


In the energy functional, there are two competing 
Casimir invariants, namely C 2 °^^\{<p)) and C2°^^\{4>))- 
By minimizing the energy functional, the following four 
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phases are obtained [IM3: 

'ferromagnetic phase : ( 0 ) = ( 0 )// = ‘ 1 , 0 , 0 , 0 , 0 ) 
for Cj < 0 and < 0 ; 

cyclic phase :((/>) = * (i, 0 , 0 , 

for c'l < 0 and C 2 > 0 ; 

uniaxial nematic phase : (^) = (</>)z = *( 0 , 0 , 1 , 0 , 0 ) 
for Cj > 0 ; 

biaxial nematic phase : {4>) = * (-^, 0 , 0 , 0 , 

. for c[ > 0 , 

( 221 ) 

where the order parameters are normalized such that 
1(0)1 = 1. We note that the uniaxial nematic and bi¬ 
axial nematic phases are energetically degenerate at the 
mean-field level. 

While the ferromagnetic phase and the uniaxial ne¬ 
matic phase are described by /r-SB with fin and /iq, re¬ 
spectively, the cyclic phase and the biaxial nematic phase 
are not. However, both the cyclic phase and the biaxial 
nematic phase are inert states. Moreover, they are both 
described by linear combinations of the highest-weight, 
zero-weight, and lowest-weight states with simple ratios 
between the coefficients, 1 and 

B. Spin-1 color superconductor 


Under the spin rotation, Ac,z transforms in the vector 
representation of ^(^(S): 

^ Ru'Ac^i, for R s SO{3). (223) 

Also, there is the f7(l)-symmetry associated with the 
baryon-number conservation which acts on the order pa¬ 
rameter A as 

A e*‘^A for 0 G R. (224) 

Based on the above discussions, the order parameter of 
the spin-1 color superconducting phase is given by a 3 x 3 
complex matrix 

A = {Ac,j|c = 1,2,3,Z = x,?/,z}, (225) 

and the symmetry group G is G = U{3) x SO{3) which 
consists of three symmetries, the color 5'f7(3)-symmetry, 
the spin S'0(3)-symmetry, and the f7(l)-symmetry as¬ 
sociated with the baryon number conservation, respec¬ 
tively. Combining Eqs. (I222p . (12231) . and (I224L the order 
parameter A transforms under G as 

At^t/*A*i?fort7GC/(3),i?GS'Cl(3). (226) 

We note that the system has a combined symmetry of 
U{3) and SO{3), resulting in the two Casimir invariants 
in the energy functional. The Cartan generators of the 
Lie algebra of G consist of three generators; two genera¬ 
tors, A 3 and As, of su(3) and one generator, Sz, of so(3). 
They are defined as 


Next, we consider the example of spin -1 color super¬ 
conductors. As we will see below, the four ground states 
are obtained from the minimization of the energy func¬ 
tional and three of them are inert while one of them is 
not. Color superconducting phases are the superconduct¬ 
ing phases in which Cooper pairs formed by quarks are 
condensed • Since each quark field ^ has three in¬ 
ternal degrees of freedom, flavor /, spin s, and color c, 
the resulting Cooper pair has these three internal de¬ 
grees of freedom. In the spin-1 color superconducting 
phases, quarks form a Cooper pair in a single flavor, a 
spin 5'0(3)-triplet, and a color S'C/(3)-antitriplet chan¬ 
nel The single flavor and the spin S'0(3)-triplet 

imply that the Cooper pair does not have an internal 
degree of freedom in the flavor but has the spin 1 , re¬ 
spectively. The color S'f7(3)-antitriplet implies that the 
Cooper pair has three color charge, anti-red, anti-blue, 
and anti-green. Let Ac,; be the field of the Cooper pair 
with color c (= 1, 2, 3) and the spin direction parallel to 
the I (= X, y, z)-axis, respectively, where colors 1 , 2 , and 
3 denote the color charge anti-red, anti-blue, and anti¬ 
green, respectively. The term “anti” implies that Ac,/ 
transforms in the conjugate representation of the three- 
dimensional representation of SU{3): 

Ac,; ^ 

when g !->• [/cc' 9 /,s {U € SU{3). ( 222 ) 


A 3 


1 0 0 \ , / 1 0 0 \ 

0 -1 0 , As = ^ P 1 0 , (227) 

oooy v3po-2y 

0 f 0 \ 

-f 0 0 . (228) 

0 0 0 / 


We note that from Eq. (12261) the actions of A 3 , As, and 
Sz on the order parameter commute because the SU (3)- 
group and its generators act on the order parameter A 
from the left, while the S'0(3)-group and its generators 
act on it from the right. The mean-field energy functional 
of spin -1 color superconduct ors can be constructed from 
the Hilbert-Schmit norm a/TL/AAT) of the matrix A 
and the Casimir invariants of SU{3) and SO{3) as 

U(A) = -cTr(AA^) +co [Tr(AA^)]% ciCr^^^(A) 
+C2G2‘’^^\a). (229) 

Here, C' 2 ^^^^(A) and C' 2 °^^^(A) are the Casimir invari- 
ants of SU{3) and 5’0(3) defined by 


8 


C'2^“'^HA)=^[Tr(AAaAt)]^ 

0—1 

(230) 

C2^“(^^(A)=^[Tr(A5,At)]^ 

(231) 


a—1 
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where {Aa}®^! andJ^a}o=i are the set of the Gell-Mann 
matrices of su(3) and the set of the generators of 
the vector representation of so(3) defined by {Sa)bc = 
respectively. In the energy functional, there are 
two competing Casimir invariants, namely C' 2 ”*'^^(A) and 
^so( 3 )(^)^ These Casimir invariants are related to the 
quartic invariants Tr(AAlAAl) and Tr [A‘A(A*A)1] 
used in Ref. [l^ as 

= 2Tr(AAl'AA^) - - [Tr(AA1■)]^ 

O 

(232) 

C 2 °^^\a) = Tr(AA^AA^) - Tr [A*A(A*A)t] . 

(233) 

From the analysis of the quartic invariants in Ref. [ 13 , 
they satisfy the following inequalities 

0<C2^“("^(A)<^, (234) 

0<Cf^^^(A)<l (235) 

for a 3 X 3 matrix normalized as Tr(AAl) = 1. By min¬ 
imizing the energy functional, the following four phases 
are obtained M- 


1 


1 i 0 


A phase : A^ = -—*10 
2 ' 0 00 ^ 

for 2ci + C 2 < 0 and C 2 < 0, (236) 

/o 0 o' 

polar phase : Ap =000 

Vo 0 1 , 

for Cl < 0 and C 2 > 0, (237) 

1 0 0 

color — spin — locked phase : Ac = —)= (010 

0 0 1 


V3 


for Cl > 0 and ci -|- C 2 > 0 , 

( ei *ei 0 
—*ei Cl 0 
0 0 £2 

for 2ci -I- C 2 > 0 and ci -I- C 2 < 0, 


(238) 


(239) 


where ei and £2 are the constant values defined as 


€1 = 


2ci 


2 V 4ci -I- C2 


£2 = 


2ci -I- C 2 
4ci -I- C 2 


(240) 


The order parameters are normalized such that 
Tr(AAt) = 1. 

Let us analyze these phases from the viewpoint of /r- 
SB, the inert state, and the Casimir invariants. Since 
the symmetry group of this system is no longer a sim¬ 
ple Lie group, we generalize the concept of /r-SB for 
ODLRO to the case in which the symmetry group of 
the system is a semisimple Lie group. Since a semisim¬ 
ple Lie algebra can be decomposed into the direct sum of 


the one-dimensional commutative Lie algebras and sim¬ 
ple Lie algebras, we can assign the Casimir invariants to 
each simple Lie algebra. When the order parameter is a 
simultaneous eigenstate of all of the Cartan generators 
and each Casimir invariant is minimized or maximized 
for the state, we refer to such a symmetry breaking as 
/i-symmetry breaking. Among the four ground states, 
the A phase, the polar phase, and the color-spin-locked 
phase are inert states, while the e phase is not. The order 
parameter of the A phase is a simultaneous eigenstate of 
A 3 , As, and S^: 


AsA^ — A^, AsA^ — , AaSz — Aa- (241) 

We note from Eq. (I226p that the generator Sz of so(3) 
acts on the order parameter from the right. In the A 
phase, both maximized: 


CT^^\Aa) = I CI°^^\Aa) = I. (242) 

Therefore, the A phase is described by ^-SB. The order 
parameter of the polar phase is a simultaneous eigenstate 
of A 3 , As, and Sz, 


AsAp = —^Ap, AsAp = 0, ApSz = 0, (243) 
V3 

and 02 “^^^A) is maximized , whereas C 2 ‘’^^^(A) is min¬ 
imized: 

C2“^^\Ap) = C' 2 '’^^V^p) = 0. (244) 


Therefore, the polar phase is described by fi-SB. The 
color-spin-locked phase is an inert state but is not /r- 
SB. However, we can see from Eq. (I238p that the ratios 
between the components are all simple numbers similarly 
to the case of the spin-2 BECs; they are all one. In the 
color-spin-locked phase, both and 

are minimized: 


Cf (^V^c) = 0, C2'^^\Ac) = 0. (245) 


The £ phase is not an inert state. This phase is an inter¬ 
mediate phase between the A phase and the polar phase. 
In the limit ci/c 2 —>■ —1/2 (ci/c 2 —>■ 0), it coincides 
with the A phase (the polar phase). In the e phase, 
the Casimir invariants takes intermediate values between 
their minimum and maximum: 


^su(3) 


^so(3) 



8ci(2ci -I- C 2 ) 
(4ci -I- C 2 )^ 
2ci y 
4ci -I- C2 y 


(246) 

(247) 


In spin-I color superconductors, a non-inert state 
emerges as a consequence of the competition between dif¬ 
ferent Casimir invariants. 
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VIII. CONCLUSION 


In conclusion, we have proposed a Lie-algebraic ap¬ 
proach to systematically finding mean fields of quantum 
many-body systems on the basis of the dynamical sym¬ 
metry. The mean fields of /r-symmetry breaking is de¬ 
rived through the minimization of the energy functional 
constructed from the Casimir invariants. We have intro¬ 
duced a concept of /r-symmetry breaking as a phase that 
is characterized by a weight vector in the representation 
of the Lie algebra. For ^-SB, the quadratic part of an 
effective Lagrangian of NG modes is block-diagonalized 
as in Eq. (I59|) in terms of the Cartan canonical form. In 
^-SB there appear three types of NG modes as listed in 
Table n Also, homotopy groups of topological excita¬ 
tions are calculated systematically for /i-SB as summa¬ 
rized in Eqs. m and (IMl) . The textures of NG modes 
and topological excitations are described in terms of the 
generalized magnetization Sa- By applying /x-SB to a 
C/(iV)-symmetric system, we have demonstrated that fi- 
SB involves a large class of symmetry broken phases as 
listed in Table ITTll 

In Sec. I VIII we have seen from the examples of spin-2 
BEGs and spin-1 color superconductors that many of the 
ground states obtained by the minimization of the en¬ 
ergy functional are inert, despite the fact that there is a 
competition between different Gasimir invariants. More¬ 
over, these states are described by linear combinations 
of weight vectors with simple ratios between the coeffi¬ 
cients. The physics behind this fact is yet to be fully 
understood and merits further study. 
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Appendix: Proof of Theorem [3l on the homotopy 
groups in p-SB 


In this Appendix, we prove Theorem[3]for /x-SB on the 
basis of the theory of an integral lattice and a co-root 
lattice (jllldll. 

Let r and R+ be the rank and the set of positive roots 
of a compact Lie group G. We define an integral lattice 
La for a compact Lie group G and a lattice Lji{S) for a 


subset S of R+ as follows [4l|, : 


Lg := {f e K 
LniS) := Span^ 


exp(27rxi7t) = e}, 
a G S 


(a, a) 


(A.l) 

(A.2) 


where Span^A denotes a vector space spanned by ele¬ 
ments of X with integer coefficients: 




1 y^nkXk 

[ fe=l 

xfe e A, rzfc e Z, zx' e N > 


Under a addition of r-dimensional vectors, Lq and Lr^S) 
of R+ become Abelian groups for any subset S. It is 
known that Lr^S) is an Abelian subgroup of Lq for any 
subset S of i?+ [i^ 113. 

The proof of Theorem [3] proceeds in four steps. 

First, we prove the following lemma on the general 
formulas of homotopy groups. 


Lemma 1 Let z* : 7r„(i7) — >• 7 r„(G) be the induced ho¬ 
momorphism of the inclusion map i : H ^ G. The homo¬ 
topy groups of the homogeneous space TTi{G/H) {i = 1,2) 
are calculated as follows: 

(1) For a Lie group G and its connected subgroup H, 

7ri{G/H) = Coker{it : 7ri(iL) ^ 7ri(G)}, (A.4) 

where Goker{/ : A —K} for a homomorphism / : A —>■ 
Y is defined as 


Coker{/ : X ^ Y} ■.= Y/lm{f :X ^Y}. (A.5) 

Let [a] (a € 7 ri(G)) and 0 € [0,27r] be a representative 
element of the coset space Goker i\ and the azimuth an¬ 
gle around the vortex. The texture O{0) of the vortex 
associated with [a] is given by 

0{e) = a{9) o Go, (A.6) 

where Oq is the value of the order parameter at 0 = 0 

and g o Oo {g G G) is the action of g on Oq- 

(2) For a compact Lie group G and its subgroup H, 

Tr2{G/H) = KerK : 7ri(i7) ^ 7ri(G)}. (A.7) 

Let tr and ae (0 < 0 < tt) be an element of Ker zj and 
a continuous deformation from ag=o = cr to (T 6 »= 7 r = e 
(trivial loop). The texture 0{9,(j)) of the point defect 
associated with a is given by 

0{9,4)) = ag{4>)oOQ, (A.8) 

where 0 G [0, tt] and (p G [0, 27r] are the three-dimensional 
polar coordinates and Oq is the value of the order param¬ 
eter a.t 9 = <j) = 0. 


Proof of Lemma [7] 

We first prove the formula 


7Tn{G/H) 

Coker z* 


= Ker z:_i 


(A.9) 








21 


by using the following homotopy exact sequence [58j |: 


^ 712 (G) 772 (G/iJ) 7ri(iJ) 7ri(G) 

^7ri(G/iI)^^o(ir)^^o(G), (A.IO) 

where p* : 7r„(G) -)> 7r„(G/iJ) and 5* : 7r„(G/i?) -)> 
7r„_i(G/-ff) are the induced homomorphism of the pro¬ 
jection map p : G ^ G/H and the boundary map 
d : G/H H. 

By using the homomorphism theorem and the exact 
sequence, we obtain 


nn{G/H) 
Ker d* 


= Im 


Im d* = Ker 


(A.ll) 

(A.12) 


Ker d* = Im p* = _ Coker j* . (A.13) 

” ” Kerp* ” ^ ' 


Combining Eqs. (IA.11I) - (IA.13I) . we obtain Eq. (IA.9I) . For 
a connected subgroup H, we have ttq{H) = 0 and hence 
Ker 7 q = 0. Thus, we obtain Eq. (IA.4I) . Let O and 
[a] (a G 7ri(G)) be an element of 7ri(G/iL) and the el¬ 
ement of Coker i\ that corresponds to O in Eq. (IA.4I) . 
Let Go is a value of the order parameter. The projection 
p : G ^ G/H coincides with an action on Oq: 


p{9)=9°Oo. (A.14) 


In fact, any element h in iL is mapped by p to Oq) p(/i) = 
Oq because an element of H acts on Oq trivially. Since 
pI : 7ri(G) —7> ni{G/H) is the induced homomorphism 
of the projection map p : G ^ G/H, the element O of 
'Ki{G/H) corresponding to the element of [ a] ( a G 7ri(G)) 
is obtained by the projection p as follows 


0{9) = p[a{0)] = a{0) o Oq. (A.15) 

Thus, we obtain Eq. (IA.6I) . Since 772(G) = 0 for a com¬ 
pact Lie group G , we obtain Coker *2=0 and hence 
Eq. (IA.7I) . Let cr and ag (0 < 0 < tt) be the element 
of Ker and the path of a continuous deformation from 
ag=o = cr to the trivial loop crg^^n- = e. We can obtain the 
element O of 'K 2 {G/H) that corresponds to cr in Eq. (IA.7I) 
as follows: 


O{0,(j)) = ag{(j)) o Oq, (A.16) 

where 0 G [0, tt] and (p G [0, 27r] are the three-dimensional 
polar coordinates and Oq is the value of the order pa¬ 
rameter at 0 = (/) = 0. Thus, we obtain Eq. (IA.8I) . which 
completes the proof of Lemma [TJ 

Second, we prove Eq. (IMl) on the second homotopy 
group. Let Lh be an integral lattice of H. We define the 
co-root lattice L_r,g for G by 

2a 

(a, a) 

It is known that the co-root lattice q is a subgroup of 
the integral lattice Lq and that the coset space L p /L b n 
is equivalent to the first homotopy group tti (G) |41| . l42l | : 

-KiiG) = Lq/L ji,G- (A.18) 



a G R-i 


(A.17) 


Also, we define L^.h as the co-root lattice of the sub¬ 
group H. For /i-SB, Lh is a subgroup of Lq and 
Lh,h is a subgroup of Lfi,G- Writing an element of 
771 (i7) = Lh/Lh,h StS t + LH,H{t G Lh), the inclusion 
map 7* : tti{H) —>• 77i(G) satisfies 

il{t + Lhh) — t + Lh,g- (A.19) 

From Eq. (IA.7I) . we have 
tt 2 {G/H) = Kerji* : 7ri(i7) —7ri(G)} 

= {* + Lr^h]^ G LH,i*i{t + Lr^h) G Lr^g} 


= {t -I- LR^H\t G Lh, t G Lfl.c} 

= {Lh r\ Lr^g)/Lr^h- (A.20) 

For p-SB with pin, Lr^q and Lr h can be rewritten from 
Theorem [T] as 

Lfl.G = Lr{R+), (A.21) 

Lr,h = Lr{R+\Rh)- (A.22) 

Therefore, we obtain the third row of Eq. (1^51) . Since 
Lr{R+) is a subgroup of Lq and all of the Cartan genera¬ 
tors are not broken ones for DLRO, we have LhL\Lr^g = 
Lr,G = Lr{R+). Therefore, we obtain the first row of 
Eq. (IMl) . For p-SB with po, Lr^g and Lr^h can be 
rewritten from Theorem [T] as 


Lr,G = Lr{R+), 

^{Lr{R+\R^) 
- \lr{R^\RO^) 


(A.23) 

for DLRO and po! 
for ODLRO and po, 

(A.24) 


and Lh is a subgroup of Lr g- Thus, we obtain the sec¬ 
ond and fourth rows of Eq. (j93p . The texture of the point 
defect in Eq. (IM)l is obtained by using Eq. (IA.8p . Let cr 
be an element of Ker i* with co-root vector 2a/{a, a). 
Here, a represents the loop on H dehned by 


<j{(j)) = exp 



for (j) G [0, 277]. 


(A.25) 


Since a G Ker i\, there exists a continuous deformation 
from a to the trivial loop. In fact, 


(Te{4>) = e 




iO-! 




for [0, fp) G [0, tt] X [0, 277] 


(A.26) 


describes a continuous deformation from ag=o = a to the 
trivial loop crg^^ = e: 


cre=o{(t>) = efo = a{p). 


(A.27) 


/ i\ —in ) ^ k id) / \ in , ^ \ 

og=T^{p) = o (“>“)e e e 
= e e = e. 


i^7 


(A.28) 


Here, we use 


e 


— in 






(A.29) 
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From Eq. (IA.8|) . the element O' £ tt 2 {G/H) correspond¬ 
ing to <j can be obtained by acting (Js{4>) on Oq: 


O'{9, (p) = e o O, 


2^7 


= o Oq, 


(A.30) 


where we use the fact that Ha is an unbroken generator. 
Through the continuous deformation 

O'^{0,4>) = e*“^ O'(0,0) for 0 < u < 1, (A.31) 

O'^ is deformed from O'^^^q = O' to 0„=i = O. Thus, 
Eq. dMl) is obtained. 

Third, we show Coker = 0 for DLRO. Erom 
Eqs. (IA.4I) and (IA.19|) . we have 

Im il = {iU-t + LR,H)\t^LH} 

= {t + Ln^clt G Lh} 

= Span^jAff U Aij.d/Afi.G, (A.32) 

and hence we obatin from Eq. (IA.18|) 


from Theorem [TJ Thus, Coker can be written as 


Coker z* 


_SpanzjJ}_ 

Span^l/} n SpanzlAff U Lr^g} 

(A.36) 


Let a„ (n £ Z) be an element of Span^j/} in the nu¬ 
merator of Eq. (IA.36E Since the action of a„(0) on the 
order parameter represents the rotation of the phase 2mr 
around the vortex, Eq. (IA.6I) reduces to Eq. (1^ . 

Eor the ODLRO with fiQ, Lr is a subgroup of Lr^g 
from Theorem [TJ We obtain from Theorem |3| 


Spanz{/} n Spanz(A_f/ n Lr^g) = Span^l/} n Lr^g 

= 0, (A.37) 

Coker z* = Span^j/}. (A.38) 

Therefore, we obtain Coker z* = Z from Eq. (IA.36I) . For 
the ODLRO with fiR, Coker zj is a subgroup of Z. To 
show that Coker zj' is a finite group, it is sufficient to 
show 


nl £ Spanx(L/z U Lr^g) for 3rz £ Z. (A.39) 
Let r be the rank of the Lie algebra g and let C = 

be the Cartan ma¬ 
trix of g, a set of simple roots of g, and the Dynkin index 
of fiR, respectively. The generator I can be written as 


Coker i\ 


_ ^ _ 

Span^lA/f U Lr^g} 


{A.33) 


For a /i-SB with DLRO, no Cartan generators are broken 
ones. Therefore, we obtain 


Lg 
Coker z* 


Lh = Lr^g, 


_ ^ _ 

Span^jAc U Lg} 


= 0 . 


(A.34) 

(A.35) 


Finally, we calculate Coker i\ for a /z-SB with ODLRO. 
For a /z-SB with ODLRO, the Cartan subalgebla gq of 
G is included in f) except for one direction Span^l/} — Z 




E 

1=1 


1 ^ ^ 


2ff 


0 ) 


(a(A, alA) 
(A.40) 


Since |/z_fzp/ — H^„ and are the elements of 

Lr U Lr g and the coefficients are all rational numbers, 
the greatest common divisor of the denominator of the 
coefficients satisfies Eq. (IA.39|) . The positive integer I 
in Eq. (El) is determined to be the minimum number 
that satisfies Eq. (IA.39I) . which completes the proof of 
Theorem |3| 
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